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INTRODUCTION 


The optimum design of structures subject to a single harmonic load has been 
considered by Icerman [1] and Plaut [2]. Icerman [1] studied a rod acted upon by a 
single harmonic tip load, and he minimized the amplitude of the tip's steady state 
displacement subject to fixed volume. His derivation of the necessary optimality 
criterion required a steady state extension to the principle of minimum potential 
energy. This new extremum principle can be used to establish the global sufficiency 
of the optimality condition, provided sufficiently restrictive conditions are 
placed upon the admissible designs. 

Plaut [2] generalized Icerman' s work by minimizing the amplitude of the steady 
state deflection at any specified location of the structure. His problem allowed 
for several harmonic loads provided all were driven at the same frequency. He 
derived the necessary optimality condition by first extending the principle of 
mutual stationary potential energy to the steady state. He did not address whether 
the optimality condition was also sufficient for the optimal design. 

This paper initially addresses a simplified version of Icerman' s problem. The 
nature of the restrictive conditions that must be placed on the design space in 
order to ensure an analytic optimum are discussed in detail. Icerman 's problem is 
then extended to include multiple forcing functions with different driving fre- 
quencies. And the conditions that now must be placed upon the design space to en- 
sure an analytic optimum are again discussed. An important finding is that all 
solutions to the optimality condition (analytic stationary design) are local optima, 
but the global optimum may well be non-analytic. 

The more general problem of distributing the fixed mass of a linear elastic 
structure subject to general periodic loads in order to minimize some measure of the 
steady state deflection is also considered. This response is explicitly expressed 
in terms of Green's functional and the abstract operators defining the structure. 

The optimality criterion is derived by differentiating the response with respect to 
the design parameters [3]. The theory is applicable to finite element as well as 
distributed parameter models. 
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AN ELEMENTARY FORCED VIBRATION 


Consider a rod of length L whose cross-sectional area S is piecewise constant. 
In the simplest case S - for 0 < x < L/2 and S = S2 for L/2 < x< L (Fig. 1). 
The design ratio S 1 /S 2 is denoted by R. The tip load is F cos a)j)t. The steady 
state response can be represented by U(x) cos wpt, where |U(L)| is the amplitude of 
the tip response. Figures 2 and 3 show, respectively, the natural frequency of the 
design and tip amplitude as a function of the design ratio R. For designs with R= 
R*, the design and driving frequencies coincide and hence the resonance condition 
shown, Icerman’s [1] optimality condition is sufficient provided only designs R >R 
are considered. 



Figure 1. Two-Segment Rod Subject To A Single Driving 
Frequence 
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FREQUENCY AND TIP-DEFLECTION AMPLITUDE VS DESIGN 


Rauwcr A 


Non-Dimensional Frequency vs. Design. The non-dimensional 
frequency A is defined as uiL/p/E where P is mass per unit length 
and E is Young’s modulus. 


Fig. 3. Tip Deflection vs. Design. R m i n is the analytic local optimum. 

Depending upon the value of U(0), R m i n may or may not be globally 
the optimum. 






TWO ELEMENT ROD SUBJECTED TO TWO DRIVING FREQUENCIES 

The rod shown in Figure 1 is now subjected to a tip load consisting of two 
frequencies - F=F 1 cos + F^ cos2 co^t. While the addition of the second driving 
frequency does not affect the period of the response, it does add a minor complica- 
tion to the calculation of the maximum response. The response can be expressed in the 
form U^(x) costop t + U 2^ x ) cos 2u>pt. The maximum tip response U max is obtained by 
setting x=L and maximizing the magnitude of the response over one period. 

* * 

There are two resonant frequencies and, in general, two designs R^ and R£ for 
which the fundamental frequency equals one of the resonant frequencies. A typical 
response is depicted in Figure A. Clearly, there are up to three, local extrema; 
all are minima and one is not analytic. Analogous to the simpler case of only one 
forcing frequency, one or both of the resonant designs may not exist. 


IUWI/L 



fl tsi/sz) 


Figure A. Typical Maximum Response as a Function of the Design Parameter. 
R^ and R^ are analytic local optimum designs. 
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OPTIMIZATION USING GREEN'S FUNCTIONAL 


Consider the following boundary value problem 

(T* ET + otM)u = f in Q 

Byu = g in 30^ 

B y ET u = h in 30^ 

* 

where T, T adjoint operators 

E(S) ... Linear stiffness operator 

M(S) ... Linear mass operator 

a A scalar 

•k 

Y, Y ... Trace operators mapping functions defined in Q onto functions 
defined on 30^ and respectively. 

>v 

B,B ... Boundary operators 

.Complementary subsets of SQ 
S Design variable(s) 

The following integration by parts formula, due to Oden and Reddy [4], is postulated 
(Tu, ETv)^ - (u,T ETv)^ 

= (YU,B*VETv) an - (Byu,y K ETv) 9 ^ 

Also E,M satisfy 

(u,Ev) n = (v,Eu) fi 

(u,Mv) = (v,Mu) for every admissible u and v. 
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STRUCTURAL DYNAMICS 


The equations governing linear structural dynamics can be expressed as 
ft 

T ET u + Mu = f in Q 

Byu = g on 
ft ft 

By ET u = h on 

If f,g, and h all have the same periodicity, 

N 

f = E f^Cx) C0S n °t) t 
N 

g = E g^x) cos n w D t 
N 

h - E h (x) cos n u) n t 

i n D 

The solution u (x,t) satisfies 

N 

u ■ E U (x) cos n w_.t 
i " D 

The Fourier coefficients satisfy 


L = (T*ET - n 2 
n 

w 2 M) U 
D n 

- f 

n 

in 

a 

B Y U n 

= 

8 n 

on 

asjj 

* 

B y ET U 
1 n 

= 

h 

n 

on 

3C2 2 


It can be shown that the solution is 

* 


U n “ “n'Vn + ( ®n' Y "VsOj + <h n*^n ) 8B J 


where G is the Green's function corresponding to the operator L . 
n n 
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SOLUTION TO VARIATIONAL EQUATIONS 


Define 

6U = U (x,S+6S) - U (x,S) 
n n 9 n 

Then, by using Oden and Reddy's integration by parts formula and taking the first 
variation of the differential equation specifying U^, it can be shown that 

6U = - (6E TU , TG ) 0 + nV (6MU ,G ) 
n nnSt D n n a 

Consequently, 

N 2 2 

6u(x,t) = I [ n w"' ( SMU , G ) 0 - ( 6ETU , TG ) 0 ] cos n w t 
j »D n n n nnw -D 

To determine whether any particular design which causes a stationary response 
at fixed x and t also minimizes the response can be obtained by variational calcu- 
lus. The second variation <S 2 u clearly involves the terms 6^E, 6U n ,6^M. However, 

<$U is determined above, and 
n 

6G = - (TG 6E TG ) n + n 2 w 2 (6MG ,G ) n 
n n n Q D n n Q 

2 

Consequently, 6 u is completely determined by the variations in the operators E 
and M. 
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APPLICATIONS TO A ROD 


A rod with varying cross-sectional area A(x) is fixed at x=0 and acted upon by 
the horizontal periodic force 

E F cos n u^t 
n D 

at the tip x=L. The steady state equations are 

- (E A U 1 ) ' - c‘ 2 E A u> 2 n 2 U = F 6(L-x) 
on oD nn 

U (0) = E A U ' (L)=0 
n on 

-2 

where E q is Young’s modulus and c = p/E. Clearly 

U (x) = F G (L,x) 
n n n 

Upon making the following identification for the operators: 

T = - T* = d/dx , E=E q A=S(x) , M=pA = c" 2 S(x) 
the solution for 6 U (L) becomes 


6U (L) = — 

n F 

n 




2 2 -2 2 
n wT c U } 6S dx 
D n J 


Now, let T be selected so that 


u(L,T) = U 


max 


By imposing the fixed mass constraint 
L 

/ <5S dx = 0, 

o 


the optimality condition for minimizing U 
tained, i.e. 


max 


over all admissible designs S is ob- 


N 

Z 

n=l 


cos n co^T 


F 

n 




-2 . 2 
c U 

n 


) = A 


For rods consisting of piecewise constant cross-sections S^. in the domain < 

x < x^, the optimality condition takes the form 


N 

Z 

n=l 


cos n (OpT 
F 

n 



x 


i-1 



2 2 -2 tt 2 . , 
n c U . Jdx 
D ni 


A. 

l 
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RESULTS FOR A TWO -SEGMENT ROD WITH ONE DRIVING FREQUENCY 

Table 1 shows results for a single driving. The non-dimensional frequency G 
is ojpL/c, the analytic optimum design stiffness ratio R is S-/S 2 , U q is the 
analytic minimum, U un £f is the amplitude of vibration for R=I and U p is the non- 
analytic minimum obtained by letting R approach zero. 

It is observed that for low frequencies, the analytic optimum is clearly supe- 
rior to any other design. For G > 0.854, the awkward non-analytic optimum provides 
a smaller amplitude for the response than does the analytic design. For still 
larger frequencies G > 1.743, even the uniform design is preferable to the analytic 
local optimum. 


TABLE 1. Comparison of Analytic Optimum to Non-Analytic Optimum 
and Uniform Designs. 


G 

R 

U ,/U 
opt umf 

U ,,/U 
opt 0 

0.1 

1.005 

0.99999 

0.010 

0.3 

1.046 

0.999 

0.093 

0.5 

1.130 

0.996 

0.278 

0.7 

1.266 

0.982 

0.604 

0.854 

1.415 

0.957 

1.000 

1.0 

1.597 

0.911 

1.625 

1.4 

2.419 

0.497 

4.066 

1.743 

3.828 

1.000 

13.657 

1.8 

4.176 

2.824 

19.734 


1274 


RESULTS FOR A TWO- SEGMENT ROD WITH TWO DRIVING FREQUENCIES 


Table 2 shows numerical results for two driving frequencies. The symbols Q 
and R are defined the same as in the previous example. Here, U q , U ^ an< * 
are, respectively, the maximum absolute value of the displacement over one periSd 
for a local analytical optimal design, uniform design, and a rod for which S- tends 
to zero. The forcing function is F^ cos oj^t + F^ cos 2o) D t. The non-dimensional 
frequency is fixed at 0.6 for this example. Rote that there are two analytic 
local optimal designs R- and R^- In the third and fourth columns in Table 2 are 
comparisons based upon the design R^ while the sixth and seventh columns are com- 
parisons based upon R^. 


The global optimum is obtained by searching the fourth and seventh columns. 
Since the entries in Column 4 always exceed the corresponding entry in Column 7, it 
follows that R 2 always provides a better solution than R- . Further, for F^ > F^> 

R 2 is the global optimum but for < F 2 f is t ^ ie Sl°Dal optimum. 


TABLE 2. Comparison of Both Analytic Optimal Designs with Non-Analytic 
Optimal Designs with Non-Analytic Optimum and Uniform Designs. 
All data is for = 0.6 


’ /f 
V 2 

R 1 

V 

U . 
opt 

R 2 

U opt 

u . 

opt 



U unif 

U 

0 


U unif 

u — 

0 

1/9 

0.168 

0.81 

2.02 

1.89 

0.80 

1.99 

2/8 

0.19A 

0.99 

1.88 

1.85 

0.81 

1.5A 

3/7 

0.213 

1.15 

1.72 

1.79 

0.83 

1.2A 

A/6 

0.23A 

1.29 

1.56 

1.7A 

0.85 

1.03 

5/5 

0.250 

1.A2 

1.A1 

1.68 

0.87 

0.87 

6/A 

0.270 

1.53 

1.27 

1.61 

0.89 

0.7A 

7/3 

0.288 

1.63 

1.13 

1.5A 

0.92 

0.6A 

8/2 

0.313 

1.69 

1.00 

1.A5 

0.9A 

0.55 

9/1 

0.3A7 

1.69 

0.85 

1.3A 

0.97 

0.A8 
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ADDITIONAL RESULTS FOR TWO DRIVING FREQUENCIES 


The identical problem as the preceding one is considered. Here, however, re- 
sults are presented for fixed force amplitudes where F^=F 2 * The driving frequency 
ftp is varied in Table 3. 

Of the analytic solutions, it is clear that R^ is superior to R^ for values of 
ftp closer to 1.0, but R^ is better for the smaller driving frequencies. However R^ 
never provides the optimum global design. For forcing frequencies ftp < 0.6, R 2 
is the global optimum; otherwise the global optimum is not analytic. 


TABLE 3. Comparison of Both Analytic Optimal Designs with 


Non-Analytic Optimum 
for F 1 /F 2 = 1. 


s 

R 1 

v 

U . 
opt 



U unif 

U 

0 

0.2 

0.025 

16.6 

1.10 

0.4 

0.105 

3.93 

1.27 

0.6 

0.251 

1.42 

1.40 

0.8 

0.482 

0.14 

1.76 

1.0 

0.824 

0.96 

2.36 


and Uniform Designs. All data is 


R 2 

U o P t 

U . 
opt 


U unif 

U 

0 

1.05 

0.999 

0.07 

1.24 

0.985 

0.30 

1.68 

0.87 

0.87 

2.67 

0.18 

2.31 

5.09 

2.84 

7.00 
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TWO SEGMENT ROD SUBJECT TO THREE FREQUENCIES 


The two segment rod is now subjected to the tip force cos u)^t + F„ cos 
2 Wpt + F^ cos 3u)pt. There are now three resonant designs and three local optimal 
designs in addition to the non-analytic optimal design as R-K). Table 4a contains 
results for the lowest driving frequency 0^ = 0.3. In all cases R„ provides the 
global optimum. However, if 0^ = 0.6 (Table 4b), then R^ is the global optimum. 


TABLE 4. Comparison of Three Analytic Local Optimum Designs with Non-Analytic 
Local Optimum Design, = 0.3 for (4a) and 0^ = 0.6 for (4b), 


V F 1 

F 3 /F l 

R 1 

U . 
opt 

U 

o 

R 2 

U . 
opt 

“U 

o 

R 3 

U . 
opt 

U 

o 

0.50 

0.50 

0.065 

0.82 

0.157 

0.40 

1.22 

0.18 

0.75 

0.75 

0.061 

0.89 

0.158 

0.52 

1.25 

0.21 

1.00 

1.00 

0.058 

0.94 

0.158 

0.54 

1.27 

0.34 





(4a) 





0.50 

0.50 

0.286 

1.07 

0.825 

0.67 

3.12 

1.38 

0.75 

0.75 

0.267 

1.15 

0.830 

0.76 

3.24 

1.72 

1.00 

1.00 

0.254 

1.20 

0.838 

0.95 

3.32 

2.02 


(4b) 
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